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Similar to the optimal-doped, weak-ferromagnetic (WFM induced by canted antiferromagnetism, 
Tcurie — 131 K) and superconducting (T c = 56 K) RuSr2GdCu208, the underdoped RuSr2EuCu20s 
(Tcurie = 133 K, T c = 36 K) also exhibited a spontaneous vortex state (SVS) between 16 K and 
36 K. The low field (±20 G) superconducting hysteresis loop indicates a weak and narrow Meissner 
state region of average lower critical field B°r(T) = B^ e (0)[l - (T/T svs ) 2 ], with B% e (Q) = 7 G 
and Tsvs = 16 K. The vortex melting transition (Tmelting = 21 K) below T c obtained from the 
broad resistivity drop and the onset of diamagnetic signal indicates a vortex liquid region due to the 
coexistence and interplay between superconductivity and WFM order. No visible jump in specific 
heat was observed near T c for Eu- and Gd-compound. This is not surprising, since the electronic 
specific heat is easily overshadowed by the large phonon and weak-ferromagnetic contributions. 
Furthermore, a broad resistivity transition due to low vortex melting temperature would also lead 
to a correspondingly reduced height of any specific heat jump. Finally, with the baseline from the 
nonmagnetic Eu-compound, specific heat data analysis confirms the magnetic entropy associated 
with antiferromagnetic ordering of Gd 3+ (J = S = 7/2) at 2.5 K to be close to NaIc ln8 as expected. 

PACS numbers: 74.72.-h, 74.25.Ha 



I. INTRODUCTION 

Anomalous physical properties have been observed 
recently in the weak-ferromagnetic (WFM induced 
by canted antiferromagnetism) and high-T c supercon- 
ducting RuS^RC^Os system (Ru-1212 with R = Sm, 
Eu, Gd, and Y) having a tetragonal TlBa2CaCu207-type 

structure !- 2 -3,4, 5, 6, 7, 8, 9 , 10, 11, 12, 13, 14, 15, 16, 17. 18, 19, 20, 21. 22, 23, 24 , 25 

Possible superconductivity was also reported in Ca- 
substituted WFM compounds RuCa2RCu2 0s (R = 
Pr-Gd)i 49 ' 50 ' 51 The weak-ferromagnetism in these 
strongly-correlated electron systems originates from the 
long range order of Ru moments in the RuOg octahedra 
due to a strong Rxi-Ad X y^ yz ^ zx -0-2p x ^ y _ z hybridization 
with a Curie temperature Tcurie ~ 131 K. A G-type 
antiferromagnetic order probably occurs with Ru 5+ 
moment /i canted along the tetragonal basal plane, 
even through the small net spontaneous magnetic 
moment fi s <§C /x(Ru 5+ ) is too small to be detected in 
neutron diffraction^! 5 -^ : 10 ! 22 The Ru valence of 4+ and 
5+ was determined from x-ray absorption near edge 
measurements ! 23 ' 52 

With its quasi-two-dimensional Cu02 bi-layers sep- 
arated by a rare earth layer in the Ru-1212 struc- 
ture, RuSr2GdCu208 has the highest resistivity-onset 
temperature T c ~ 60 K among different Ru-1212 



compoundsiii 2 ^' 5 -' 3 ^ A broad resistivity transition width 
AT C = T c (onset) - T c (zero) = T c - T meltmg ~ 15-20 K 
is most likely a consequence of coexistence and interplay 
between superconductivity and WFM order. The dia- 
magnetic signal is observed only near Tmelting instead 
of T c , and a reasonable large Meissner signal can be 
detected only in zero-field-cooled (ZFC) mode. 47 Lower 

and Sm-compound7respectivelyj ' la JNo superconductiv- 
ity can be detected in RuSr 2 RCu 2 8 (R = Pr, Nd)^ 
while a superconducting RuSr 2 YCu20s phase is stable 
only under the high pressure i 21 : 26 



Interest of the current work stimulates from a recent 
report of spontaneous vortex state (SVS) between 30 K 
and 56 K in RuSr2GdCu 2 0s- 47 However, the compound 
undergoes a low temperature antiferromagnetic order- 
ing arising from Gd 3+ at 2.5 K. To avoid this compli- 
cation, isostructural RuS^EuC^Os with nonmagnetic- 
Eu 3+ ions was chosen as a prototype material in this 
study to evaluate the anomalous magnetic, transport, 
calorimetric properties and d-wave nature near and be- 
low T c = 36 K. The calorimetric data were further used 
as a basis in elucidating the magnetic entropy associated 
with the Gd 3+ ordering. 
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II. EXPERIMENTAL 

Stoichiometric RuS^RCv^Os samples were synthe- 
sized by solid-state reactions. High-purity R11O2 (99.99 
%), SrC0 3 (99.9 %), R 2 3 (99.99 %) (R = Pr, Nd, Sm, 
Eu, and Gd), and CuO (99.9 %), in the nominal com- 
position ratios of Ru:Sr:R:Cu = 1: 2: 1: 2, were well 
mixed and calcined at 960°C in air for 16 hours. The 
calcined powders were then pressed into pellets and sin- 
tered in flowing N2 gas at 1015°C for 10 hours to form 
RuSr2RC>6 and CU2O precursors. This step is crucial in 
order to avoid the formation of impurity phases. The 
N2-sintered pellets were heated at 1060°C in flowing O2 
gas for 10 hours to form the Ru-1212 phase, then oxygen- 
annealed at a slightly higher 1065°C for 7 days and slowly 
furnace-cooled to room temperature with a rate of 15°C 
per hour4i 

Powder x-ray diffraction data were collected with a 
Rigaku Rotaflex 18-kW rotating-anode diffractometer us- 
ing Cu-K Q radiation. Four-probe electrical resistivity 
measurements were performed with a Linear Research 
LR-700 ac (16Hz) resistance bridge from 2 K to 300 
K. Magnetic susceptibility and magnetic hysteresis mea- 
surements from 2 K to 300 K in low applied magnetic 
fields were carried out with a Quantum Design //-metal 
shielded MPMS2 superconducting quantum interference 
device (SQUID) magnetometer. Calorimetric measure- 
ments were made from 1 K to 70 K by using a thermal- 
relaxation microcalorimeter. A mg-size sample was at- 
tached with a minute amount of grease to a sapphire 
holder to ensure good thermal coupling. The sample 
holder had a Cernox temperature sensor and a Ni-Cr al- 
loy film heater. The holder was linked thermally to a 
copper block by four Au-Cu alloy wires. The tempera- 
ture of the block could be raised in steps but held con- 
stant when a heat pulse was applied. Following each heat 
pulse, the sample temperature relaxation rate was moni- 
tored to yield a time constant r. The total heat capacity 
was calculated from the expression c = kt, where k is the 
thermal conductance of Au-Cu wires. The heat capacity 
of the holder was measured separately for addenda cor- 
rection. The molar specific heat of the sample was then 
obtained from C = (c - c a( Menda)/(m/M) with m and M 
being the sample's mass and molar mass, respectively. 



III. RESULTS AND DISCUSSION 

Figure 1 summarizes structural and superconducting 
properties, as a function of R 3+ ionic radius r (coordina- 
tion number CN = 8), of various RuS^RC^Os-a system 
(R = Pr-Y). T c decreases from a maximum value of 60 
K for optimal-doped Gd (r = 0.105 nm) to 36 K for un- 
derdoped Eu (r = 0.107 nm), and < 10 K for Sm (r = 
0.108 nm). Larger rare earth ions of Nd (0.112 nm) and 
Pr (0.113 nm) lead to a metal-insulator transition. Pow- 
der x-ray Rietveld refinement study indicates that the 
insulating phase is stabilized in the undistorted tetrag- 




FIG. 1: The variation of superconducting transition T c and 
tetragonal lattice parameters a, c with rare earth ionic radius 
R 3+ (coordination number CN = 8) for RuSr2RCu20g_i sys- 
tem (R = Pr-Y). 



onal phase (space group P4/mmm) with a larger lattice 
parameter a ~ 0.390-392 nm, which gives a reasonable 
Ru 5+ -0 bond length of d ~ 0.197 nm if the oxygen con- 
tent is slightly deficient (S > 0). On the other hand, the 
metallic phase with smaller rare earth ions can be sta- 
bilized in the full-oxygenated (5 ~ 0), distorted tetrago- 
nal phase (space group P4/mbm) with smaller a/y/2 ~ 
0.383-0.385 nm but still a reasonable Ru-0 bond length 
through RuOe octahedron rotation. 

Indeed, the powder x-ray diffraction pattern for the 
oxygen-annealed RuSr2EuCu20s-5 sample indicates sin- 
gle phase with tetragonal lattice parameters of a = 
0.5435(5) nm and c = 1.1552(9) nm. A Raman scatter- 
ing peak of 265 cm" 1 indicates that the Ai g mode sym- 
metry belong to a P4/mbm instead of P4/mmm group. 
Accordingly, with RuO@ octahedra rotation angle 9 ~ 
14° around the c-axis and oxygen parameter S ~ 0, 10 
Rietveld refinement analysis with a small residual error 
factor R = 5.31% yields a reasonable Ru-0 bond lengths 
d = (a/2y/f)(l - sin^)" 1 / 2 = 0.198 nm. It is close to 
the minimum calculated bond length d(Ru 5+ -0) of 0.197 
nmiiS 

Figure 2 shows the temperature dependence of field- 
cooled (FC) and zero-field-cooled (ZFC) volume mag- 
netic susceptibility 4:nxv a t TG for bulk and pow- 
der RuSr 2 EuCu208 samples. Weak-ferromagmagnetic 
ordering occurs at Tcurie = 133 K. Similar to 
RuS^GdC^Osr^ this Eu-compound has its electrical 
resistivity data, which are also included in Fig. 2, exhibit- 
ing a non- Fermi-liquid- like behavior above Tcurie- The 
linearly temperature-dependant values of 10.0 mf2 cm at 
300 K and 5.5 mf2 cm at 160 K give an extrapolated value 
of 2.6 mtt cm at K, yielding a ratio p(300 K)/p(0 K) of 
3.9. Below T curie, a T 2 behavior prevails. The onset of 
deviation at 36 K from such a temperature dependence 
is taken as the superconducting transition temperature 
T c . The melting temperature of superconducting vortex 
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FIG. 2: The electrical resistivity p(T) and volume magnetic 
susceptibility 4-7rxv(T) in 1-G field-cooled (FC) and zero- 
field-cooled (ZFC) modes for oxygen-annealed bulk and pow- 
der RuSr2EuCu208 samples. 
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FIG. 3: The low-field superconducting hysteresis loops M- 
B a at 2 K for RuSr 2 GdCu 2 08 and RuSr 2 EuCu 2 8 . Average 
lower critical field B c i(ave) at peak values and afe-plane B°t 
for deviation from initial linear lines are indicated by arrows. 



liquid is assigned to ^melting — 21 K, where resistivity 
reaches zero4£ The broad transition width of 15 K is the 
common feature for all reported Ru-1212 compounds. It 
indicates that the superconducting Josephson coupling 
along the tetragonal c-axis between Cu-0 bi-layers may 
be partially blocked by the magnetic dipole field ~Q dipole. 
of ordered Ru moments in the Ru-0 layer.— 

The Meissner shielding at 2 K is complete (47rxv = 
47rM/B a ~ 1.3) for ZFC bulk sample, but much reduced 
(-0.1) in the powder sample. However, in 1-G FC mode, 
no such an effect can be detected below T me i t i ng due to 
strong flux pinning. 

Low-field (±20 G) superconducting hysteresis loop at 2 
K for bulk sample RuSr 2 EuCu20s and RuSr 2 GdCu20s 
as reference are shown in Fig. 3. The initial magne- 
tization curve deviates from straight line at 2 G and 

3 G for the Eu- and Gd-compound, respectively. The 
narrow region of full Meissner effect roughly reflects the 
temperature-dependent lower critical field in the aft-plane 
B^(T). The average lower critical field B^" 3 for bulk sam- 
ple as determined from the peak of initial diamagnetic 
magnetization curves is 7 G for R = Eu and 13 G for 
R = Gd. The effect on the exact peak value due to the 
surface barrier pinning is neglected. For RuS^EuC^Os, 
B™ e decreases steadily from 7 G at 2 K to 6 G at 5 K, 

4 G at 10 K, and below 1 G at 15 K. A simple empirical 
parabolic fitting gives B™ e (T) = B^ e (0)[l - (T/T S vs) 2 }, 
with average B^ e (0) ~ 7 G and spontaneous vortex 
state temperature Tsys = 16 K. The Ginzburg-Landau 
anisotropy formula B^ e = (2B"^ + B^-J/S, then provides 
an estimated c-axis lower critical field B^) ~ 17 G and 
anisotropy parameter ~ 8.5. 

The lower field superconducting phase diagram for the 
polycrystallinc bulk sample is shown in Fig. 4. The 
average lower critical field B™ e separates the Meissner 
state and vortex state. The upper critical field B C 2 and 
vortex melting field R r neiung determinated from magne- 
toresistivity measurements are field-independent below 
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FIG. 4: The low field, low temperature superconducting phase 
diagram B a (T) of RuSr 2 EuCu 2 Os. The spontaneous vortex 
state (SVS) occurrs between Tgys = 16 K and T c = 36 K. 
Vortex lattice/glass melting temperature T melting is denned 
from temperature at which resistivity drops to zero. 



20 G. The WFM-induced internal dipole field B dlpole of 
8.8 G on the Cu02 bi-layers is estimated using extrapo- 
lated B a c r value at T = 0, (B^ e («) + B dlpo , e )/B^ e (0) 
= Tc/^svs- It further yields a small net spontaneous 
magnetic moment fi s of 0.1 hb per Ru, based on the 
relation of Bdipole ~ 2/i s /(c/2) 3 , where c/2 = 0.58 nm 
is the distance between midpoint of Cu02 bi-layers and 
two nearest-neighbor Ru moments. If the WFM struc- 
ture is indeed a G-type antiferromagnetic order with 1.5 
Hb for Ru 5+ in t 2g states canted along the tetragonal 
basal plane, the small /i s would give a canting angle of 
4° from the tetragonal c-axis and be difficult to be de- 
tected in neutron diffraction with a resolution ~ 0.1 fiB- 
The molar specific heat data up to 70 K in Fig. 5 
show a good agreement between Eu- and Gd-compounds, 
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FIG. 5: The molar specific heat of RuSr2RCu208 (R = Eu, 
Gd). Antiferromagnetic Gd 3+ ordering prevails at 2.5 K. 

except that a peak reflects the antiferromagnetic Gd 3+ 
ordering near T^y ~ 2.5 K. Consistent with previous re- 
sults for lower-T c Gd-compounds in zero applied mag- 
netic field i 15 ' 28 No visible jump in specific heat was ob- 
served near T c = 36 K. This is not surprising, since only 
the electronic component in specific heat would change 
with superconducting transition, but it is easily overshad- 
owed by the much larger phonon contribution. Specifi- 
cally, assuming a same magnitude as that observed in 
Lai.gsSro.isCuCU (AC ~ 0.33 J/mol K at T c = 37 K) 
and YBa 2 Cu 3 7 (AC - 4.6 J/mol K at T c = 92 K), 53 
an estimated AC ~ I J/mol K at T c here is only about 
1% of total specific heat, falling below the experimen- 
tal precision. In addition, the broad resistivity transition 
due to vortex melting would further points to a corre- 
spondingly reduced height of AC. 

It would be of interest to obtain information on the 
Gd 3+ ordering. To do so, delineation of various contri- 
butions to the total specific heat begins with the non- 
magnetic Eu-compound up to 7 K. In the format of C/T 
versus T 2 , the data in Fig. 6 can be well fitted by the sum 
of four terms with different temperature dependence: 

C = (3T 3 + aT 2 +jT+ 7 ^r. (1) 

The coefficient of the first term, /3 = 0.89 mJ/mol K 4 , 
can be used to derive a Debye temperature #d of the 
lattice, 

f3 = n(l2ir 4 /5)N A k/9 3 D , (2) 

where Na is Avogadro's number, k the Boltzmann con- 
stant, and the number of atoms per formula unit n = 
14. The 6d value of 312 K thus obtained supports the 
validity of the T 3 -dependence approximation in Debye 
model for the lattice specific heat below 7 K ~ 6jj/50. 
The quadratic term has two possible sources: the nodal 
line excitation for d-wave pairing symmetry and the spin 
wave excitation of WFM Ru sublattice. The fact that the 
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FIG. 6: Low temperature C/T versus T 2 of RuSr 2 EuCu 2 8 
from 1 K to 7 K. Data above 1 K can be fitted using C(T) = 
/3T 3 + aT 2 + 7T + 7//T 2 with Debye temperature 6 D = 312 
K. 



observed a value of 4.2 m J/mol K is much large than 0.1 
mJ/mol K of YBa2Cu307 could be an indication of a less 
important nodal line excitation, but an enhanced spin 
wave excitation. The linear term is considered normally 
as an electronic contribution, which is not expected to ex- 
ist in a superconductor at temperature much lower than 
T c . While the observed coefficient 7 = 7.3 mJ/mol K 2 is 
comparable to that of some cuprates, its origin remains 
to be identified. One plausible explanation is based on 
the complicated magnetic structure and mixed valence. 
Such a scenario could lead to a spin glass-like lattice, for 
which an even larger linear term in specific heat has been 
observed in another Ru compound of Ba 2 PrRu06- 54 

The last term with aT" 2 dependence is most likely the 
high-temperature tail of a Schottky anomaly. Its occur- 
rence at the relatively low temperatures suggests nuclear 
energy splittings being the cause. Such energy splittings 
occur typically for nuclei having a spin I and magnetic 
moment fi n in a hyperfine magnetic field ffhf- For the 
calorimetrical measurements under consideration, they 
are is most likely associated with the Ru nuclei, since the 
Ad magnetic moments of ordered Ru ions are spatially 
fixed, polarizing the s-electrons and producing a net spin 
at the nuclei, yielding a hyperfine field. There are two 
Ru isotopes with non-zero /i„: "Ru (fractional natural 
abundance A = 0.1276, I = 5/2, and (j, n = -0.6413) and 
101 Ru (A = 0.1706, I = 5/2, and fx n = -0.7188)^ How- 
ever, nuclear energy splittings can also be caused by the 
interaction between the qudrupole moment of a nucleus 
and the electric field gradient produced by neighboring 
atoms. The electric field gradient could be quite high in 
the layered compound. Meanwhile, Cu and Eu or 155 Gd 
(A = 14.7%) and 157 Gd (A = 15.7%) nuclei all have non- 
zero quadrupole moment. Without the full knowledge 
of magnetic hyperfine field and electric field gradient, it 
is not feasible at present to delineate the experimentally 
obtained 77 of 6.63 mJ K/mol into the two different con- 
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FIG. 7: Temperature dependence of magnetic specific 
heat and entropy (inset) associated with Gd 3+ ordering in 
RuSr 2 GdCu 2 08. 

tributions. 

By assuming that its various coefficients in Eq. (1) for 
Eu-compound remain the same for the Gd-compound, 
One can then obtain the magnetic contribution to specific 
heat associated with antiferromagnetic Gd 3+ ordering as 

C m = Cad — Ceu- (3) 

The results are shown in Fig. 7. Using the format of 
C TO /T versus T. It is of interest to note a broad shoul- 
der below T^v, a common feature seemingly prevailing in 
other similar type of compounds such as GdBa2Cu307, 
GdBa 2 Cu 4 8 and TlBa 2 GdCu 2 7 ^^^ According to 
Fishman and Liu^ it is due to spin fluctuations in the 
normally ordered state, and such fluctuations are more 
pronounced for large spins. Indeed, Gd 3+ has the largest 
spin among all R 3+ ions. The areal integral in Fig. 7, 
including that associated with the broad shoulder should 
yield the magnetic entropy, 

S m = f (C m /T)dT. (4) 



As shown in the inset, S m reaches a saturation value 
of 17.6 J/mol K around 10 K. Considering the built-in 
approximation in Eq. (4), it agrees exceptional well with 
the theoretical value of N A k ln(2J+l) = N A k ln8 = 17.2 
J/mol K for the complete ordering of Gd 3+ . 



IV. CONCLUSION 

The lower critical field with B cl (0) = 7 G and Tgys 
= 16 K indicates the existence of a spontaneous vortex 
state (SVS) between 16 K and T c of 36 K. This SVS state 
is closely related to the weak-ferromagnetic order with 
a net spontaneous magnetic moment of ~ 0.1 /is/Ru, 
which generates a weak magnetic dipole field around 8.8 
G in the Cu0 2 bi-layers. The vortex melting transition 
temperature at 21 K obtained from resistivity measure- 
ments and the onset of diamagnetic signal indicates a 
broad vortex liquid region due to the coexistence and in- 
terplay between superconductivity and WFM order. No 
visible specific heat jump was observed near T c for Eu- 
and Gd-compound, since the electronic specific heat is 
easily overshadowed by the large phonon contributions 
and the expected jump would spread over a wide range 
of temperature due to vortex melting. Finally, the mag- 
netic entropy associated with Gd 3+ antiferromagnetic or- 
dering at 2.5 K is confirmed to be close to N A k ln8 for J 
= S = 7/2. 
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